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Abstract 

We present a model which predicts the temperature of Bose-Ein- 
stein condensation in atomic alkali gases and find excellent agreement 
with recent experimental observations. A system of bosons confined 
by a harmonic oscillator potential is not characterized by a critical 
temperature in the same way as an identical system which is not con- 
fined. We discuss the problem of Bose-Einstein condensation in an 
isotropic harmonic oscillator potential analytically and numerically 
for a range of parameters of relevance to the study of low temperature 
gases of alkali metals. 
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A stunning development in the last year has been the demonstration of 
Bose-Einstein condensation (BEC) in gases of rubidium 0, lithium 0, and 
most recently sodium This experimental work has stimulated theoretical 
studies to try to understand the underlying physics of the situation [|, [5], §]. 
For many purposes the complicated magnetic traps used in the experiments 
can be approximated by harmonic oscillator potentials. There have been 
several studies of BEC in harmonic oscillator confining potentials f?], . 
The purpose of our paper is to examine the condensation of bosons in a 
harmonic oscillator potential in a detailed way, and to point out a serious 
flaw in the analyses of Refs. || [9|, [TO]; namely, unlike the situation for 
a boson gas with no external confining potential |IT |, there does not exist 



a critical temperature which signals a phase transition. However, we will 
show that there is a temperature at which the specific heat has a maximum 
which can be identified as the temperature at which BEC occurs. By using 
parameters relevant to the experiments with alkali gases |1|, |2], |3jwe find 
excellent agreement with the experiments. 

The fact that a gas of bosons (neglecting interactions) in a harmonic oscil- 
lator potential does not have a phase transition at some critical temperature 
is already apparent from the early work of 0. (This will also be shown be- 
low in a very simple way.) A similar situation occurs for a system of charged 
bosons in a homogeneous magnetic field; in three spatial dimensions BEC 
does not occur in the same way as for the case where there is no magnetic 



field [0. The same is true for bosons confined by spatial boundaries |13 



A general criterion to decide whether or not BEC occurs has been given re- 



cently by us |14|, and it is easy to show that the criterion is not met for a 
system of bosons confined by a harmonic oscillator potential. 

Since the experiments of Refs. [[[], @, || obviously are seeing BEC, a 
natural question which arises concerns the exact nature of the phenomenon. 
If BEC as found normally for the free boson gas is impossible for a system 
of bosons in a confining potential, then in what sense does BEC occur ? A 
natural criterion which has been used in systems of finite size has been to 
look at the maximum of the specific heat ]T3| . We will apply this criterion 
to the case of bosons confined by a harmonic oscillator potential. 

The simple model which we will discuss here is appropriate to neutral 
spin-0 bosons confined by a harmonic oscillator potential. For mathematical 
simplicity we will only deal with the case of an isotropic potential in this 
paper. ( The more general analysis of anisotropic potentials, along with 
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mathematical details relevant to the analysis presented here will be given 
elsewhere |16| .) We will assume that the system can be described by a grand 
canonical ensemble. The grand potential is defined by 

q=-J2ln(l-zexp(-(3E N )) (1) 

N 

where (3 = (kT)" 1 , E N are the energy levels, and z = is the fugacity 
in terms of the chemical potential /i. It proves convenient to expand the 
logarithm in ([1]) to obtain 

oo z n 

g = E-E ex p(-^^). (2) 

n=l n N 

For an isotropic harmonic oscillator characterized by an angular frequency 
uj the energy levels are given by E% = (k + 3/2)hu> with multiplicity (k + 

l)(k + 2)/2 where k — 0, 1,2, The sum over N in (|2|) may be performed 

to obtain 

q = £ -(1 - e-)3 (3) 



where we have defined the dimensionless variable x = hu/(kT). The number 



of particles is given hy N = [3- 1 [^-) , which becomes 



N 



7> 

oo ( ,nP(p—3/2huj) 



n=l 



In order that the number of particles remain positive, it is necessary for 
H < 3/2hu). (More generally, we require \i < E where E Q is the lowest energy 
level.) Normally the critical temperature for BEC is the temperature at 
which n = E Q , which for the isotropic harmonic oscillator reads fi = 3/2tkv. 
It is now easy to see that BEC cannot occur in the same way for bosons 
confined in the harmonic oscillator potential as it does for bosons in free 
space. In the case of the free boson gas in free space with no confining 
potential, as the temperature is lowered the chemical potential /i increases 
from negative values towards the value 0. (This is in agreement with the 
general result [i = Eq quoted above since the lowest energy level is zero for 
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the free boson gas.) The value of the temperature at which \x = defines 
a critical temperature T c determined in terms of the particle density. At 
temperatures lower than T c , \x remains frozen at the value \x = 0, and the 
number of particles found in excited states is bounded. If the total number 
of particles exceeds this bound then the only possibility is for the excess 
particles to be found in the ground state, giving rise to BEC. This standard 



scenario is described in []IT] in some detail. 

In the case of (f|) it is easy to see that the number of particles satisfies 
the following simple bound (assuming w ^ 0) 

e /3(n-3/2hu) 

N > l _ e /3(Ai-3/2M ' ( 5 ) 

As \x — > 3/2hu, which should signify the onset of BEC, this expression shows 
that N increases without bound, unlike the free boson gas where as \i — > 
TV has an upper bound. To phrase this another way, no matter how low T 
becomes (so long as it remains non-zero) it is always possible to solve (|J) for 
\i with n < 3/2Tiuj for any finite value of iV no matter how large N is. If N, 
as given in (^) had an upper bound, as for the case of the free boson gas in 
three spatial dimensions, this would no longer be the case. There is therefore 
a qualitative as well as a quantitative difference between the free boson gas 
and a system of bosons confined in a harmonic oscillator potential. 

We now come to our analysis of ([$]) and (|J) and our fundamental difference 
with earlier treatments || |], |10J- We have treated, and will continue to treat, 
the energy spectrum for bosons confined in the harmonic oscillator potential 
as discrete. Refs. || |9|, [UJ make the physically plausible but mathematically 



unsound assumption that if x = fun/ (kT) << 1, then it is possible to replace 
the sums which arise by integrals. It is precisely this procedure which leads 
to the erroneous conclusion that a critical temperature exists at which (in 
our units) fi = 3/2hu. Obviously if the original sum (^) has no upper bound 
as \i — > 3/2fiu, any approximation which leads to an upper bound is suspect. 
If the behaviour for small x is desired, a safer approach is to take the exact 
results expressed in terms of sums, and to try to obtain expansions valid for 
small x. This is also useful for numerical purposes since the original sums 
(0) and dU) are not very rapidly convergent if x is small. For sodium ||, if 
we take T ~ 2 x 10~ 6 K, and take u to be given by the geometric mean of the 
frequencies of 235 Hz, 410 Hz, and 745 Hz (i.e. u;/27r=4i6 Hz), then x is of 



4 



order 10~ 2 . Thus the behaviour of the thermodynamic quantities for small 
x is of interest. 

In addition to the dimensionless variable x defined earlier, it is convenient 
to define a new dimensionless variable e related to the chemical potential \i 
by jj, = hu>(3/2 — e). If we assume that x « 1 and e << 1, then the leading 
behaviour for q in addition to N is well accounted for by approximating 
e~ nex (l — e~ nx )~ 3 ~ (nx)~ 3 which is obtained by a simple Taylor expansion. 
(For the anisotropic oscillator this also shows that identifying uo with the 
geometric mean of the frequencies is justified for the leading order term.) 
This approximation can be criticized on the grounds that no matter how 
small we take x, n eventually becomes large enough so that nx » 1. A 
more satisfactory procedure from the mathematical standpoint is to convert 
the sums defining q and N into contour integrals, and to deform the contour 
in an appropriate way for obtaining at least an asymptotic expansion of the 
desired function. The details of this procedure are somewhat involved and 
will be given elsewhere [16]. 



If x « 1 and e << 1, then it can be shown that 

N ~ C(3)x- 3 + (ex)- 1 + §C(2)^- 2 - — - 6C(2)aT 2 (6) 

is a good approximation. (((z) denotes the Riemann ^-function.) It is possi- 
ble to obtain an expansion for q to the same level of accuracy; however, the 
expression involves more complicated functions than that for iV and will not 
be given here. It can be noted that the leading term in (|6|) agrees with that 
found using the naive procedure mentioned in the previous paragraph. It is 
possible to obtain an expression for the specific heat C valid for x « 1 and 
e « 1 : 

k~ l C ~ 12C(4)x- 3 + yC(3)x- 2 . (7) 

The internal energy is given by 

U~hu [3C(4)a;- 4 + -((3) X - 3 + -((2) X - 2 + -(ex)- 1 - 3eC(3)aT 3 j . (8) 

These expressions agree to a good approximation with numerical results 
found by analyzing the exact sums when x << 1 and e << 1 hold. In 
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particular, the result for N in (^) is useful for determining the chemical 
potential in terms of the number of particles and the temperature. 

For large x (in fact if x > 1) the original sums (^j and (f|) converge quite 
rapidly and are useful for numerical computation. For x » 1 the specific 
heat becomes 

k~ x C ~ 3x 2 e~ x . (9) 

The specific heat therefore vanishes as x — > oo corresponding to T — > 0. It is 
also possible to obtain analytic approximations valid for x « 1 but without 
assuming e << 1. These expressions all involve more complicated functions 
than those we have presented here and will not be given. 

Unfortunately it is difficult to find a reliable analytic approximation for 
the specific heat which allows us to study whether it has a maximum, and if 
it does, at what temperature it occurs. Instead we have studied the specific 
heat numerically by doing the relevant sums for a range of temperatures. It 
is necessary to solve (f|) for the chemical potential in terms of N. When 
this is done for the case N = 5 x 10 5 [Q, it is found that over a very small 
range of x, corresponding to a very small range of temperatures, e changes 
from large values of the order of 100 to small values of the order 10~ 3 . This 
is a good indicator that the specific heat has an unusual behaviour in this 
temperature range. We have shown our results for the specific heat over the 
range of interest in Fig. 1. 

The ordinate gives the specific heat in units of the Boltzmann constant k, 
and the abscissa is x = %u/{kT) with toj'l'n = 416 Hz. It is readily apparent 
from this figure that the specific heat undergoes a very rapid change over a 
very small range of temperatures. The maximum in the specific heat occurs 
for x ~ 0.0136 which corresponds to a temperature of T ~ 1.47 x 10~ 6 K. 
This is in remarkable agreement with the value of 2 x 10~ 6 K in Ref. [[|. 
It supports our proposal that the BEC observed in sodium gas at least is 
associated with the maximum in the specific heat rather than with a phase 
transition as occurs for the free boson gas. 

The specific heat is perfectly continuous and smooth at its maximum. 
We have shown the behaviour of the specific heat over a small range of x 
where it changes most rapidly in Fig. 2. Our results contrast markedly with 
the discontinuous behaviour found in Refs. || |9], [TIJ which was found by 
approximating sums with integrals. This further supports our claim that 
such an approximation leads to unreliable conclusions. 
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Figure 1: The specific heat in units of k as a function of x = hco/(kT). 
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Figure 2: The specific heat in units of k as a function of 100a; where x 
hu/{kT). 
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Although we have concentrated on sodium, we can easily extend the anal- 
ysis to rubidium and lithium || . We have worked out the maximum in the 
specific heat for N = 2 x 10 5 and N = 2 x 10 4 , and find that the maximum 
occurs at x ~ 0.0185 and x ~ 0.0408 respectively. (The precise graphs in 
these two cases are qualitatively the same as those in Figs. 1-2.) Using the 
oscillator frequencies given in Ref. || for lithium, and taking N = 2 x 10 5 we 
find the temperature at which the specific heat has a maximum is T ~ 380 
nK in good agreement with the range of 100-400 nK for the experiment. For 
rubidium with iV = 2 x 10 4 we find T ~ 71 nK if we use the oscillator fre- 
quencies of Ref. [|I| . If we use the frequencies given in Ref. j| for the strong 
trap we find T ~ 124 nK, again in close agreement with experiment. Given 
the simplicity of our model, the predicted condensation temperatures are in 
remarkable agreement with experiment. 

In conclusion, we have studied the condensation of spin-0 bosons confined 
by an isotropic harmonic oscillator potential. We have shown that although 
there is no phase transition such as that occurring in the free unconfined 
boson gas, it is possible to identify a temperature at which BEC occurs by 
looking at the maximum in the specific heat. Using parameters relevant to 
the experiments which have been performed on alkali gases, we find that our 
predicted temperatures are in very close agreement with experiment. Al- 
though the gases in the experiments are dilute, it would be of interest to 
include weak interactions among the atoms, since interactions have proven 
important in other contexts j|, ||, ||. Our results seem to indicate that the 
transition temperature is not affected much by interactions. In any case, 
it would not be prudent to use the free boson gas as the zeroth order ap- 
proximation in any perturbative treatment of interactions. Finally, it is of 
extreme interest to compare our theoretical predictions for the behaviour of 
the specific heat with experiment in a much more substantial way. 
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